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SUMMARY 


A review  is  presented  of  recent  results  of  ultrasonic 
investigations  of  rare  earth-iron  compounds  and  alloys.  The  discussion 
is  divided  into  three  sections.  The  first  describes  results  of  sound 
velocity  measurements  of  single  crystal  Tb  3Dy^7Fe2  as  a 
function  of  magnetic  field  strength  and  orientation  with  respect 
to  the  crystallographic  axes.  In  the  second  part  recent  work  on 
polycrystalline  and  amorphous  samples,  including  evidence  for  simul- 
taneous propagation  of  two  different  types  of  magnetoelastic  waves 
is  presented.  A model  developed  to  describe  this  phenomenon  is 
outlined  as  well. 

Finally  a transmission  experiment  in  a single  crystal  (Tb^ 3Dy^ 7Fe2) 
is  described.  Resonant-like  transmission  is  shown  to  take  place 
when  an  applied  magnetic  field  is  so  oriented  as  to  produce  almost 
equal  velocities  for  the  and  niodes. 
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PREFACE 


The  work  reported  here  represents  a portion  of  the  research 
carried  out  by  the  Magnetism  and  Materials  Group,  Solid  State  Branch 
Research  and  Technology  Department,  to  determine  the  elastic  proper- 
ties of  highly  magnetostr ictive  rare  earth-iron  alloys.  The  work 
was  supported  by  the  Office  of  Naval  Research  (Task  No.  RROll-08-01) 
The  report  summarizes  results  on  single  crystal,  polycrystalline 
and  amorphous  alloys. 

J.  R.  DIXON 
By  direction 
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INTRODUCTION 

Magnetically  ordered  solids  generally  undergo  a change  in  length 
relative  to  their  length  in  the  demagnetized  condition.  This  rela- 
tive length  change,  or  magnetostriction  is  an  intrinsic  property 
of  the  magnet  and  varies  in  size  from  nearly  one  percent  in  ferro- 
magnetic Rare  Earth  metals  at  low  temperatures  to  a few  parts  in 
a million  in  the  ferromagnetic  transition  elements  Co,  Fe  and 
Ni.  [1]  Magnetostriction  results  from  the  coupling  of  the  magnetic 
moments  to  the  lattice  distortions.  This  same  interaction  is  capable 
of  modifying  the  elastic  properties  of  magnetostr ictive  materials. 
Besides  the  changes  in  modulus  brought  on  by  domain-wall  movements 
such  as  the  AE  effect  12] , there  are  often  measurable  changes  in 
the  elastic  properties  in  magnetically  saturated  materials.  The 
latter  effects  are  most  clearly  observable  in  single-crystal  speci- 
mens. Changes  in  elastic  constants  are  most  pronounced  in  highly 
magnetostr ictive  magnets.  In  Ni  for  example,  changes  in  elastic 
constants  of  the  order  of  10“3  have  been  observed  [3];  in  rare  earth 
metals  these  changes  can  be  several  percent  [4] . 


1.  Clark,  A.  E.,  in  Handbook  of  the  Physics  and  Chemistry  of  the 
Rare  Earths,  K.  Gschneidner  and  H.  Eyring,  eds.  North  Holland 
(Amsterdam) , to  be  published. 

2.  Bozorth,  R.  M. , Ferromagnetism,  Van  Nostrand  1951. 

3.  Alers,  G.  A.,  Neighbours,  J.  R. , and  Sato,  H. , J.  Phys,  Chem. 
Solids  9,  21  (1958) . 

4.  Moran,  T.  J.,  and  Luthi,  B. , J.  Phys.  Chem.  Solids  31,  1735 


(1970) . 
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Fig.  1 Magnetostriction  constant  determined  by  the 

x-ray  method. 


For  fer romagnets , an  estimate  for  the  change  in  elastic  constant 
can  be  made  from  the  formula 

6C/C  = -X^C/K.  (1) 


Here  C is  a modulus,  X the  magnetostriction  and  K the  magnetic 
anisotropy.  For  Ni,  C =10l2erg/cm3 , K •/'lO^  erg/cm^,  X^IO"^  and 
Equation  (1)  gives  6C/C  >/'10“3.  That  the  changes  in  modulus  are  so 
dependent  on  magnetic  anisotropy  simply  reflects  the  fact  that  the 
slowing  of  a sound  wave  is  accomplished  by  a rotation  of  the  moment. 

Recently,  attempts  have  been  made  to  produce  materials  which 
would  have  the  high  magnetostriction  characteristic  of  the  rare  earths 
but  have  this  property  available  at  room  temperature.  Success  was 
achieved  [1]  with  the  cubic  rare  earth-iron  compounds  RFe2»  where 
R stands  for  a rare  earth  element.  In  TbFe2  example, 

Xiii  = 2.6  X 10”3  at  room  temperature.  Since  this  is  two  orders 
of  magnitude  larger  than  the  magnetostriction  in  Ni,  we  would  expect 
modulus  changes  of  several  percent  in  these  compounds  were  it  not 
for  their  characteristically  large  magnetic  anisotropy.  The  latter 
difficulty  can  be  overcome  with  little  expense  in  magnetostriction, 

, however,  by  forming  pseudobinary  compounds  Ri-x^'x  with  a 

suitable  choice  of  x.  In  Tb.3Dy^7Fe2,  for  instance  K is  less  than 
10®  erg/cm3  while  X is  still  large  (1.6  x 10”^).  The  reason  this 
tailoring  of  X and  K works  is  that  these  two  quantities  depend  almost 
j exclusively  on  the  rare  earth  constituent.  Using  these  two  values, 

the  formula  gives  6C/C=-1  at  room  temperature.  This  was  more  than 
ample  motivation  for  an  ultrasonic  investigation  of  the  elastic 
constants  of  this  type  of  pseudobinary  compound. 
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One  other  characteristic  property  of  these  compounds  is  of  special 
interest  here.  That  is  that  the  magnetostriction  is  highly 
anisotropic;  ^iii>>^100*  This  property  is  displayed  in  Figure  1 
in  which  the  x-ray  distortion  is  given  as  a function  of  x in 
Tbi_xDy3jFe2 . [5]  For  samples  with  x<.7,[lll]  is  the  easy  axis  and 
the  distortion  (Xm)  is  huge.  For  x>.7,  the  samples  have  [100]  easy 
and  the  distortion  (X^oo^  disappears. 

We  expect  little  anisotropy  in  the  pure  elastic  constants 
(a  = 2 C44/Cii-C]^2  ~ 1*25  in  Tb  3Dy  •7Fe2)  . The  anisotropy  in  X then 
translates  into  anisotropy  in  the  magnetoelastic  coupling: 


In  succeeding  sections  of  this  parer  we  (a)  review  our  recent 
ultrasonic  velocity  measurements  on  single-crystal  Tb^3Dy  7Fe2  from 
which  the  elcStic  constants  and  b2  were  determined;  (b)  discuss  recent 
results  of  velocity  measurements  in  highly  magnetostr ictive  poly- 
crystalline and  amorphous  rare  earth-iron  alloys.  Finally,  we 
describe  a measurement  on  single-crystal  Tb^3Dy^7Fe2  of  a resonant 
coupling  of  the  C44  mode  to  the  Ci]^-C]^2  inode  which  is  a direct 
consequence  of  the  anisotropy  in  magnetostriction  and  the  fact  that 
6C/C  > a-1. 

SINGLE  CRYSTAL  RESULTS:  and 


The  velocity  of  transverse  sound  waves  was  measured  in  a single 
crystal  of  Tb^3Dy^7Fe2.  [6]  The  direction  of  propagation  was  [110]. 
All  our  measurements  were  done  at  room  temperature,  using  frequencies 
near  10  MHz.  Figure  2 shows  the  results  of  applying  a magnetic  field 
of  up  to  25  kOe  in  various  directions  for  both  [001]  (C44)  and  [iTO] 

(C11-C12)  polarizations  of  the  sound.  Note  that  the  velocity  of 
[001]  waves  has  been  reduced  by  thirty  percent  upon  rotating  the 
field  in  the  (110)  plane  from  perpendicular  to  the  polarization  to 
parallel  to  it.  This  corresponds  to  more  than  a 50  percent  reduction 
in  C44  as  we  anticipated  from  the  estimate  given  above.  The  velocity 
of  the  Cii-C]^2  waves  was  increased  by  a percent  or  so  as  the  magnetic 
field  was  rotated  from  perpendicular  to  the  polarization  to  parallel 
to  it.  From  the  form  of  the  magnetoelastic  energy  in  terms  of  the 
magnetization  components  M^i  and  the  strains  eij  for  cubic  systems, 
ie.  we  see  that  b2  couples^components  of  the  moment  to  C^^-type 


*^1  2 ^2 
^ (Mf-l/3)e..  + Z.  M,M^e, 

me  ..2  1 1 ' 11  „2  i/j 


M‘ 


M 


111] 


5.  Clark,  A.,  Cullen,  J.,  McMasters,  0.,  and  Callen 
Conf.  Proc.  No.  29,  (A.I.P.  New  York)  p.  129. 

6.  Rinaldi,  S.,  Cullen,  J.  R.,  and  Blessing,  G.  V. , 
61A,  465  (1977). 


(2) 


E. , A.I.P. 

Phys.  Lett. 
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Fig.  2 Sound  velocity  vs.  magnetic  field:  polarizations  along 
[001]  (circles)  and  along  [llOl  (crosses).  For  each 
polarization,  data  for  two  orientations  of  field  are  shown. 


(Gij  i/j)  strains  while  is  appropriate  for  C]^i-C]^2 
strains.  Thus  the  large  changes  in  C44  and  relatively  small'^ones  in 
confirm  the  results  of  earlier  work  that  b2>>b2* 

From  Equation  (2)  and  the  quasi-static  approximation  [7] , which 
is  good  whenever  the  sonic  frequencies  are  much  less  than  the  ferro- 
magnetic resonance  frequency,  it  is  possible  to  calculate  the  change 
in  the  elastic  energy  in  terms  of  the  susceptibility  tensor  and  the 
components  of  the  moment  and  magnetic  field.  For  H //  [001]  the 
change  in  is  predicted  to  be  (H°  is  the  total  internal  field,) 

C44  = C°4  - b2/(H+H°)Mg  (3) 

i.e.  the  sum  of  the  anisotropy  field  and  4irM3  in  the  case  of 
conducting  materials  like  these. 

Here  C°44  is  the  pure  elastic  constant.  The  change  in  C44  for 
H//11101  is  zero;  thus  C44 (H//[ 110] ) = C44  and  by  plotting 
(C44  - 044)“^  vs.  H,  |b2l  was  determined  from  the  slope  of  the  , 

resulting  straight-line  portion  (see  Figure  3)  to  be  2.3  x 10^  erg/cm"^ 
in  excellent  agreement  with  earlier  measurements.  [1] 


7.  See  for  example,  Le  Craw,  R.  C.,  and  Comstock,  R.  L. , in 

Physical  Acoustics  ed.  W.  Mason  Vol  III  B page  127  (Academic 
Press  1965) . 
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Fig.  3 The  reciprocal  of  the  difference  in  €44,  (044-044)”^  plotted 
vs.  magnetic  field.  b2  is  determined  from  the  slope. 


Further  verification  of  the  magnetostr ictive  origin  was  obtained 
from  a plot  of  (Cmax-C) / (Cmax-Cmin)  versus  the  angle  of  the  magnetic 
field  with  respect  to  the  crystalline  axes  (Figure  4) . The  high- 
field  limit  of  the  theory  [8]  based  on  (2)  is  the  solid  line  in  this 
figure.  We  see  that  the  over-all  agreement  is  quite  good,  better 
at  10  kOe  than  at  6 kOe.  This  is  understandable  considering  that 
the  theoretical  curve  is  a high-field  limit.  Returning  to  the 
(C?4-C44)-l  vs.  H plot,  from  the  intercept  and  (3)  we  obtained  the  . 
magnetic  anisotropy.  Using  800  Oe  for  Mg  we  found  K=  4 x 10°  erg/cm'^ 
about  eight  times  greater  than  the  value  obtained  from  more 
conventional  methods.  [9] 

We  mentioned  that  the  velocity  of  the  C11-C12  type  mode  actually 
increased  slightly  when  according  to  (2) , it  should  have  decreased 
(see  Figure  5) . This  is  evidence  for  a non-linear  term  in  the 
magnetoelastic  coupling.  Since  the  shift  in  Cii-C]^2  is  independent 
of  applied  field,  the  non-linear  coupling  [10]  is  probably  the 
"morphic"  type,  [11] 


E' 

me 


F 


MiMj 


e,  , e 
kl  mn 


(4) 


8.  Simon,  G. , Naturf,  Z.,  13A  84  (1958). 


9.  Williams,  C.,  and  Koon,  N.,  Conference  on  the  Rare  Earths  and 
Actinides,  J.  Phys.  (to  be  published) . 

10.  Bonsall,  L.  , and  Melcher,  R.  , Phys.  Rev.  "lA,  1128  (1976); 
Dohm,  V.,  and  Fulde,  P. , Z.  Phys.  B21 , 368  (1975). 

11.  Mason,  W. , Phys.  Rev.  715  (1951). 
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where  F is  a tensor  of  rank  six.  We  will  return  to  this  question 
of  non-linearity  in  the  final  section. 


(001)  1110)  llTO) 

11101  IITOI  IK1I  11101 


Fig.  4 Angular  variation  of  the  relative  change:  (Cmax'C) 

in  elastic  constant  (C44)  in  the  (001) , (110)  and  (110)  planes. 
Circles:  6 kOe  data.  Crosses:  10  kOe.  The  full  lines 

are  a result  of  linear  magnetoelastic  theory. 


• ••• 


(ITO)  1110)  (001) 

11101  )Doi)  nloi  11101 


Fig.  5 Angular  variation  of  the  relative  change :_  (Cmax"^)  / ^^max~*"min^ 
in  elastic  constant  (Cn  - C22)  in  the  (110)  , (110)  an(3  (001) 
planes.  All  data  taken  at  10  kOe.  Dashed  lines  correspond 
to  the  morphic  effect;  full  lines  correspond  to  linear 
magnetoelastic  theory. 
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POLYCRYSTALLINE  AND  AMORPHOUS  MATERIALS;  ELASTIC  ISOTROPY 

Changes  in  moduli  in  magnetostr ictive  polycrystalline  materials 
are  very  often  observed  and  usually  attributed  [2]  to  the  movement 
of  domains.  Changes  in  Young's  modulus  were  reported  earlier  at 
low  fields  [12]  in  rare  earth-iron  alloys.  Modulus  changes  in  a 
series  of  polycrystalline  RFe2  compounds  were  measured  at  10  MHZ 
in  fields  up  to  25  kOe.  [131  The  authors  noted  that  in  TbFe2  satura- 
tion was  not  achieved  even  at  25  kOe.  The  fact  that  such  large 
changes  were  observed  at  ultrasonic  frequencies,  too  high  for  domain 
walls  to  follow,  led  us  to  believe  that  rotation  of  the  magnetization 
by  the  strain  field  was  the  cause  of  the  effects.  We  therefore 
made  [14]  a study  of  the  sound  velocity  as  a function  of  the  magnitude 
and  orientation  of  a magnetic  field  in  a series  of  highly  magnetostr ic- 
tive polycrystals  with  low  magnetic  anisotropy.  Included  in  the 
study,  besides  Tb.3Dy^7Fe2  were  Sm^ 33Dy ^ 2^2^®2 ' S"’.7^°.3^®2 
amorphous  TbFe2  (a-TbFe2) • Figure  6 shows  the  pattern  of  echoes 
we  observed  from  a-TbFe2  taken  with  the  field  at  various  angles  6 
with  respect  to  the  direction  of  polarization  of  the  transducer. 

For  0 = 0°  or  90°  we  observed  a single  train  of  echoes,  though  those 
at  0°  were  delayed  relative  to  those  at  90°.  At  intermediate  0, 
however,  we  observed  a superposition  of  the  0°  and  90°  patterns 
indicating  that  two  shear  waves  are  propagating  simultaneously,  with 
different  velocity,  in  these  configurations.  It  is  clear  from 
Figure  6 that  the  height  of  the  shorter-time  pulse  (amplitude  of 
the  faster  wave)  decreases  as  0 decreases  from  90°,  while  the  delayed- 
pulse  height  (slow-wave  amplitude)  increases.  We  checked  our 
interpretation  of  the  double  echo  pattern  by  making  a through  trans- 
mission experiment  with  the  transmitter  and  receiver  shear 
polarizations  mutually  orthogonal  (etI^r)  ^nd  rotating  the  magnetic 
field  in  the  plane  perpendicular  to  the  propagation.  Using  a large 
pulse  width  to  allow  superposition  of  the  two  waves,  the  expected 
sin  0 cos  0 dependence  of  the  amplitude  of  the  detected  signal  v;as 
verified.  We  obtained  a qualitative  description  of  the  two  waves 
simultaneously  present  by  writing  down  a coupling  between  strain 
and  magnetization  in  a way  analagous  to  Equation  (2) , but  taking 
into  account  the  over-all  isotropy  of  amorphous  and  polycrystalline 
matter.  For  isotropic  materials  there  can  be  only  one  coupling 
constant  b,  since  the  magnetoelastic  energy  E must  be  rotationally 
invariant.  Thus 


=me  = *=  «ijl 


(5) 


12.  Savage,  H. , Clark,  A.,  and  Powers,  J.,  IEEE  Trans.  Mag. 
MAG-11,  1355. 

13.  Klimker,  H. , Rosen,  M. , Dariel,  M.  P.,  and  Atzmony,  U. , Phys. 
Rev.  BIO,  2968  (1974). 


14.  Blessing,  G. , Cullen,  J.,  and  Rinaldi,  S.,  to  be  published, 
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ECHO  ECHO  ECHO 


observed  echo  pattern  for  shear  waves  propagating  perpendi- 
cular to  the  applied  field  H in  amorphous  TbFe2-  The 
changes  in  the  amplitudes  of  the  peaks  as  the  angle  0 between 
the  direction  of  H and  the  direction  of  polarization  is 
varied  indicate  the  presence  of  two  normal  modes  with 
polarization  parallel  and  perpendicular  to  H.  The  positions 
of  the  peaks  do  not  change,  indicating  that  there  is  no 
variation  in  the  velocity  of  the  normal  modes. 
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The  Mj[  are  the  components  of  the  average  magnetization  normalized 
to  one  at  saturation.  Suppose  the  field  and  thus  the  average  moment 
to  be  along  the  z axis,  and  a sound  wave  is  propagating  along  x, 
polarized  perpendicular  to  x,  but  otherwise  arbitrary.  The  only 
non-zero  strain  components  are  e^y  and  Gxz*  order  one; 

Mx  and  My  are  of  the  same  order  as  the  strain.  From  (5)  only 
is  coupled  in  first  order  to  changes  in  the  magnetization: 


E = b M Me 
me  z X xz 


(6) 


M^  can  be  obtained  from 


Mx  = X h/M^ 


(7) 


dM. 


where  x is  the  perpendicular  susceptibility  (3]^)  and  h^  is  the 


magnetoelastic  field,  defined  by 


(8) 


Thus 


_ 9 _9  9 9 

E = - b X M e /M 
me  ^ z xz'  s 


(9) 


The  total  increment  in  elastic  energy  in  the  presence  of  the  sound 
wave  is  then 


+ ) 
xy' 


(10) 


Uq  is  the  elastic  shear  modulus.  The  factor  of  1/2  in  front  of  the 
second  term  in  (10)  comes  from  the  contribution  from  the  change  in 
magnetic  energy,  which  is  -1/2  E^jg.  From  (10)  we  see  that  the 
degeneracy  of  the  shear  modes  is  broken.  The  normal  modes  are  waves, 
linearly  polarized  along  the  z and  y axes,  i.e.,  parallel  and  perpen- 
dicular to  the  field.  These  conclusions  moreover  do  not  depend  for 
their  validity  on  the  simplified  model  just  presented.  Because  of 
elastic  isotopy  in  the  plane  perpendicular  to  the  propagation  only 
the  magnetic  field  picks  out  a special  direction,  and  we  expect  the 
normal  modes  to  be  polarized  either  parallel  or  perpendicular  to 
the  field.  This  principle  applies  to  any  elastically  isotropic 
system  in  the  presence  of  a magnetic  field;  the  size  of  the  splitting 
is  proportional  to  the  square  of  the  magnetoelastic  coupling.  Since 
b is  close  to  b2/2,  which,  as  we  mentioned  at  the  outset  is  two  orders 
of  magnitude  larger  than  the  coupling  in  any  other  room-temperature 
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magnet^  it  is  not  surprising  that  this  effect  has  never  been  reported 
before  in  polycrystals.  The  same  effect  of  course  can  take  place 
in  single  crystals  if  the  propagation  is  for  example  parallel  to 
a llOOl  direction  in  a cubic  crystal.  Again,  elastic  isotropy  of 
the  (100)  plane  can  be  broken  by  a field  in  this  plane.  This  effect, 
magnetoacoustic  birefringence,  was  predicted  [15]  essentially  using 
Equation  (3) , and  has  been  observed  in  single  crystals  of  Ni  and 
Fe304  [16].  In  [16],  velocity  differences  of  the  order  of  10”^  were 
measured  using  a phase-sensitive  technique.  We  have  observed  relative 
changes  (Av/v)  of  twenty-six  percent  in  Sm  ^Ho  3Fe2  (Table  1) . 

Magnetic  field  dependences  of  both  velocities  were  separately 
measured  to  25  kOe.  Results  are  plotted  for  all  four  samples  in 
Figure  7.  The  splitting  into  two  waves  was  clearly  observable  to 
fields  of  25  kOe,  the  largest  available. 
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Fig.  7 Velocity  as  a function  of  applied  field  for  and  H//ein 

(a)  Tb.3Dy  7Fe2,  (b)  Sm.88Dy. 

(d)  a-TbFe^.  * ’ ’ ' 


The  velocities  of  the  slow  modes  could  not  be  saturated  even 
at  25  kOe;  the  velocity  of  the  fast  mode  of  a -TbFe2  is  apparently 
saturated  above  2 kOe  (Figure  7d) . The  fast  modes  of  the  polycrystals, 
especially  the  compounds  containing  Sm  show  significant  field 
dependencies  to  25  kOe.  By  plotting  the  modulus  of  the  fast  modes 
versus  1/H  we  obtained  the  pure  elastic  modulus  Mq  from  the  extra- 
polation 1/H  -►  0 (Figure  8)  . The  values  so  obtained  are  listed  in 

15.  Kittel,  C.,  Phys.  Rev.  3^,  836  (1958). 

16.  Liathi,  B.,  Appl.  Phys.  8,  107  (1966). 
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Table  1.  Then,  from  the  slopes  of  (uQ-y)"^  vs.  H plots  (Figure  9), 
we  obtained  a pair  of  coupling  parameters  h//  and  b for  the  slow 
and  fast  modes  respectively  for  each  sample.  These  are  also  given 
in  Table  1. 


Fig.  8 Modulus  (y  ) vs.  reciprocal  of  the  applied  field  for 

Tb  3Dy^7Fe2  polycrystal.  y was  determined  from  the  velocity 
of* the’ fast  mode.  From  this  plot,  yQ  was  determined. 


TABLE  1 

Elastic  and  magnetic  data  for  the  RFe2  alloys  investigated.  The 
saturation  magnetization  Mg  is  in  emu/cm^;  V is  the  sound  velocity; 
elastic  modulus  yo  is  in  10^^  erg/cm^;  and  the  effective  magneto- 
elastic coupling  constants  b and  b^^  are  in  10^  erg/cm^. 


RFe2  Alloy 

^s 

(AV) 

'V  'max 

^^o 

b 

‘-// 

Tb  3Dy_7Fe2 

780 

0.13 

4.3 

0.13 

1.3 

0.01 

S"’.88°y.l2^®2 

360 

0.21 

3.4 

0.47 

1.4 

0.11 

Sm  7H0. 3Fe2 

145 

0.26 

3.5 

0.35 

0.9 

0.15 

a-TbFe2 

390 

0.06 

2.9 

0.07 

0.5 

0.02 
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Fig.  9 Reciprocal  of  the  difference  of  the  modulus  from  the  pure 
elastic  modulus  (vi-Uq)  vs.  magnetic  field  for  Tb,3Dy  7Fe2* 
From  the  slopes  the  b's  were  determined,  a)  soft  mo5e, 
b)  stiff  mode. 


As  far  as  the  slow  1 I £)  modes  are  concerned,  the  simple  model 
outlined  above  describes  the  data  fairly  well.  Since  for  large  H, 

X =Mg/H,  Equation  (10)  predicts  for  the  change  in  modulus  <Sm//  of 
the  slow  mode  '' 


6 = - b///MgH  (11) 

Equation  (11)  correctly  gives  the  large  H behavior  of  the  slow  mode. 
It  also  predicts  that  the  modulus  changes  are  larger,  for  a given  b, 
for  samples  with  smaller  saturation  magnetization,  in  agreement  with 
the  facts;  the  largest  reduction  in  modulus  occurs  in  Sm^yHo  3Fe2f 
which  has  the  lowest  Mg  (Table  1).  (Table  1 also  shows  that’the 
b's  are  roughly  the  same  for  all  the  polycrystals.) 
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On  the  other  hand,  the  model  is  inadequate  for  a description 
of  the  field  dependence  of  the  fast  mode.  A generalization,  based 
again  on  isotropic  magnetoelastic  coupling  (Equation  (5))  but  taking 
some  account  of  the  disorder  in  the  magnetization,  predicts  that 
there  is  a change  in  modulus  for  the  stiff  mode  6y  , but  going  as 
1/h3  for  large  H.  This  generalization  is  accomplished  by  recognizing 
that  the  requirement  of  isotropy  can  still  be  met  if,  instead  of 
averages  of  the  components  of  the  magnetization,  we  couple  the  strains 
to  products  of  the  local  magnetization  components. 

Although  an  improvement,  the  above  generalization  cannot  reproduce 
the  observed  1/H  behavior  of  the  modulus  of  the  fast  modes.  The 
key  to  further  improvement  is  the  concept  of  coupling  on  a local 
scale,  i.e.  for  each  crystallite  or  region  of  coherence  in  the 
amorphous  case.  The  magnetoelastic  energy  in  a given  cubic 
crystallite  is  just  Equation  (2)  plus  surface  terms  coming  from 
couplings,  elastic  or  magnetic,  to  other  crystallites  or  regions. 

The  local  magnetoelastic  field  has  components 


M. 


(12) 


To  second  order  in  the  strain,  the  energy  in  the  crystallite  is 
lowered  by 


6E 


hiXijhj 


(13) 


The  hi  are  components  of  a real  vector  only  when  bi=b2;  in  general 
their  value  depends  on  the  set  of  axes  defining  them.  Thus  even 
if  the  magnetization  and  the  strain  were  uniform  the  hi  would  vary 
from  crystallite  to  crystallite.  The  material  has  its  total  energy 
lowered  by  the  interaction  (Equation  (13))  at  the  expense  of  inducing 
a disorder  in  the  magnetization.  Therefore  at  very  high  magnetic 
fields,  in  the  presence  of  a strain  there  are  fluctuating  components 
of  the  magnetization  in  the  plane  perpendicular  to  the  net  magnetiza- 
tion direction;  these  in  turn  will  couple  to  the  strain  components 
in  this  plane.  This  is  the  mechanism  which  leads  to  a 1/H  dependence 
of  the  modulus  of  the  stiff  mode. 

The  arguments  given  above  should  hold  true  even  if  there  are 
interactions  between  the  grains  or  local  regions,  just  as  long  as 
these  do  not  completely  inhibit  the  creation  of  magnetization 
fluctuations  via  the  magnetoelastic  interaction.  For  the  limiting 
case  of  independent  crystallites  it  is  possible  to  calculate  the 
change  in  modulus  of  either  mode  at  high  fields.  Assuming  a uniform 
distribution  of  crystallite  axes,  6E  (Equation  (13))  is  averaged 


18 


NSWC/WOL  TR  77-180 


over  the  unit  sphere  for  a set  of  fixed  strains, 
the  changes  in  shear  moduli  are 


! 

P ; 

I i . = - 

I 


8 

105 


(b2-bi) 

HM 


s 


2 


The  results  for 

(14) 


for  polarization  perpendicular  to  the  field  (stiff  mode)  and 


6y 


// 


4 ^1~^2 


^1~^2,  2.  ^ ^ 

HM_  " " HM 


T5'b^+b2 


(15) 


for  the  polarization  parallel  to  the  field  (soft  mode).  As  we 
anticipated,  the  1/H  behavior  of  the  fast  mode  depends  entirely  on 
the  anisotropy  in  the  single  crystal  magnetostriction.  In  fact, 
the  approach  to  saturation  magnetostriction  was  previously  shown  [17] 
to  behave  in  a similar  way.  As  we  mentioned  in  the  introduction, 
it  is  characteristic  of  RFe2  materials  that  Xm>>XioO  (b2>>bi)  . 


In  the  large  anisotropy  limit,  i.e.  bi  -►  0 from  (14)  and  (15) 
the  ratio  of  the  coefficients  of  the  1/H  law  for  the  two  modes  is 
predicted  to  be 


( 


2 

b2-0 


0.163 


(16) 


In  the  limit  of  magnetoelastic  isotropy,  i.e.  bi=b2,  the  prediction 
is  that  bx/b//  = 0.  The  experimental  values  listed  in  Table  1 all 
fall  between  these  two  extremes  of  the  theory.  From  the  discussion 
following  Equation  (13)  we  take  this  to  indicate  the  presence  of 
interactions  between  grains.  However  this  intergranular  interaction 
is  not  strong  enough  to  eliminate  the  magnetoelastic  anisotropy. 
Metallographic  work  shows  that  the  samples  containing  Sm,  which 
according  to  Table  1 are  closest  to  behaving  like  collections  of 
independent  crystallites,  contain  appreciable  amounts  of  intergranular 
precipitates.  Much  finer  grain  boundaries  were  observed  in  the 
Tb  7^62 

There  are  aspects  of  this  work  on  disordered  RFe2  materials  still 
to  be  worked  out,  including  a description  of  the  low  field  behavior 
of  the  moduli  and  measurements  of  their  temperature  dependencies. 

The  most  important  result  of  this  section  and  which  leads  us  into 
the  next,  is  the  breaking  of  the  elastic  isotropy  by  the 
magnetoelastic  coupling. 


17.  Lee,  E.  W. , Proc.  Phys.  Soc.  A 67^,  381  (1954). 

1 

I 
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RESONANT  TRANSMISSION  OF  SOUND  IN  A SINGLE  CRYSTAL 

Up  to  this  point  there  has  been  a certain  dichotomy  in  the 
presentation;  on  the  one  hand  we  described  single  crystal  results 
in  terms  of  changes  in  velocity  of  modes  (C44,  C11-C12)  of  fixed 
polarization  while  on  the  other  we  interpreted  polycrystalline  data 
by  means  of  modes  of  constant  velocity  whose  polarization  depended 
on  the  direction  of  the  field.  In  describing  the  single-crystal 
results  we  have  implicitly  assumed  that  the  appropriate  symmetry 


Fig.  10  Experimental  arrangement  for  the  resonance  experiment  in 
Tb.3Dy,7Fe2  single  crystal.  A shear  wave  is_generated  by 
a transmitting  transducer  with  polarization  ct//  tOOH  and 
propagates  through  the  sample  along  the  (llOl  direction. 

At  the  other  side  of_the  sample  a receiving  transducer  with 
polarization  gr//  [llOl  X gt  detects  a signal  due  to  the 
rotation  ({>  of  the  polarizations  of  the  normal  modes 
(gi,  induced  by  magnetoelastic  interaction  with  the 

magnetization  M.  The  magnetization  is  assumed  parallel 
to  the  field  H applied  in  the  (110)  plane  at  an  angle  0 
from  the  [ 001 1 . 


for  elastic  modes  was  determined  by  the  crystal.  Clearly  there  is 
some  competition  between  the  crystalline  symmetry  and  that  induced 
by  the  magnetic  field  via  the  magnetoelastic  coupling.  Generally 
the  latter  is  taken  to  be  a perturbation  of  the  former  [3,81,  so 
that  the  directions  of  polarization  of  the  transverse  modes  are  only 
slightly  tilted  by  an  angle  (|)  from  the  crystalline  axes  for  any 
orientation  0 of  the  magnetic  field  with  respect  to  the  [OOll  axis. 
(See  Figure  10).  This  tilting  of  the  polarization  is  predicted  to 
occur  only  if  both  b2  and  bi  (Equation  (1))  are  non-zero.  In  fact 
a straightforward  but  lengthy  calculation  of  the  velocities  and 
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directions  of  polarization  of  the  normal  modes  can  be  carried  out; 
[8,  18]  the  resulting  equations  for  the  velocities  are; 


2 C+C 

+ 2 

+ ((^^)^  + B^)^ 

(17) 

The 

parameters  C,  C and  B are: 

C E 

C°4/p  - b2COS^6  / (H+H°)MgP 

(18) 

C E 

: - b^sin^0  / 

(H+H°)MgP 

(19) 

B E 

b^b2sin6cos9  / (H+H°)MgP 

(20) 

where  p is  the  density. 

If  B=0  (6=0  or_90°)  there  are  two  transverse  waves  polarized 
along  [001]  and  [110]  with  velocities  and  C'^/^  respectively; 

though  either  C or  C is  reduced  (depending  on  whether  0 = 0 or  90°) 
by  the  magnetoelastic  coupling,  they  are  nonetheless  "pure"  modes 
in  that  they  are  polarized  along  crystalline  axes.  In  general  we 
only  expect  strong  effects  if  C=C'  assuming  bi  and  b2  t 0.  We  know 
from  just  looking  at  the  curves  of  Figure  (2)  that  this  will  happen 
for  some  0=6c  in  Tb_3Dy_7Fe2  because  C<C'  for  0=0  but  C>C'  for 
0 = 90°.  (0  is  measured  with  respect  to  [001]).  Furthermore,  b]^<<b2 

in  RFe2  compounds  so  the  range  of  0 around  6(,  where  strong  coupling 
can  take  place  is  small.  We  have  plotted  v versus  0 in  Figure  (11) 
using  known  values  of  b2r  Mg,  C44,  C]^]^-C]^2  taking 

b;^^  = 0.1  b2»  H=8  kOe.  To  see  what  happens  to  the  polarization 
direction  under  these  conditions,  we  calculated  ())  from 


4 

i 


taniji  2 2 Ik 

(C'-C)/2  + ((C'-C)  /4  + B^)^ 


(21)  ! 

I 

I 


4)  is  very  close  to  0 or  90°  except  in  the  immediate  vicinity  of  0 
where  it  changes  abruptly.  ° 

To  check  this  "pulling"  of  the  axes  of  polarization  off  the 
crystalline  axes  we  performed  a transmission  experiment  (Figure  (10)) 
in  which  receiving  and  transmitting  transducers  were  arranged  with 
linear  polarization  along  crystallographic  axis  and  orthogonal  to 
one  another.  A magnetic  field  of  8 kOe  was  rotated  in  the  (110) 
plane,  and  the  output  of  the  receiver  recorded.  We  expect  this 
output  to  vary  as  A sin  2^,  where  A is  another  oscillatory  function 
of  0 coming  from  interference.  Using  the  same  set  of  parameters 
listed  above,  we  have  calculated  the  expected  amplitude  as  a function 
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Fig.  11  Dispersion  in  the  velocity  of  the  shear  normal  modes  for 

[110]  propagation  in  Tb^3Dy^7Fe2  caused  by  the  interaction 
of  the  strain  with  the  magnetic  moment  in  the  (110)  plane 
at  an  angle  9 from  the  [001] . The  parameters  used  are: 
bi  = 2 X 10°  erg/cm^;  b2  = 2.3  x 10^  erg/cm°;  (C11-C32)/ 

2 = 3.9  X lO-*--^  erg/cm^;  C44  = 4.8  x lO-*-^  ezg/cm^; 

Mg  = 1000  emu/cm°;  p =9.2  g/cm°,  Hq  = 20  kOe , H = 8 kOe. 

of  0.  The  result  is  plotted  in  Figure  (12a).  The  envelope,  peaked 
around  45°,  is  due  to  sin  2(J)  while  the  oscillations  come  from  A. 

We  have  also  plotted  a modified  version  of  the  amplitude  obtained 
by  correcting  Figure  (12a)  for  depolarization  and  finite  resolution: 
Figure  (12b).  The  experimental  results  are  given  in  Figure  (12c). 
There  is  reasonably  good  agreement,  i.e.  the  theory  reproduces  the 
sharp  central  peak  and  much  of  the  structure  away  from  the  crossover 
angle.  Although  this  agreement  verifies  our  picture  of  the  rotation 
of  the  polarization,  the  theoretical  model  must  still  be  made 
consistent  with  our  earlier  discussion  of  Figure  (2)  in  which  we 
pointed  out  that  the  slight  stiffening  of  the  C11-C32  mode  was 
inconsistent  with  Equation  (2) . Including  non  linear  terms  of  the 
form  of  Equation  (4)  means  simply  replacing  bj/Ms(H+HQ)  by  itself 
plus  f]^  and  bib2/Ms  (H+Hq)  by  itself  plus  f2  in  Equation  (19)  and 
(20)  [18]  . f]^  must  be  negative  and  larger  in  magnitude  than 
b?/Mg(H+HQ)  to  explain  the  stiffening  of  C33-C]^2'  Since  the  relative 
change  in  C11-C22  is  about  2%,  f^  must  be  about  10^®  erg/cm^.  Now 
the  width  of  the  central  peak  in  the  transmission  amplitude  is  best 
fit  by  choosing  b3b2/Mg (H+Hq)  + f2  = 10^^  erg/cm^.  If  |fi|=lf2l 
must  be  the  combination  bj^b2  that  is  responsible  for  coupling  tlie 
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to  be  published. 
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Fig.  12  Comparison  between  the  theoretical  prediction  and  the 

experimental  results  of  the  resonance  experiment  described 
in  the  text  and  in  Figure  10.  The  parameters  used  for  the 
theoretical  curves  are  the  same  as  in  Figure  11. 


tv;o  transverse  modes.  At  this  stage  of  our  knowledge,  we  can  only 
regard  these  estimates  as  tentative. 

The  expectation  of  large  changes  in  elastic  properties  of  RFe2 
compounds  due  to  moment  rotations  has  been  amply  justified.  The 
measurements  described  here  have  further  shown  that  the  elastic 
isotropy  of  polycrystalline  and  amorphous  RFe2  alloys  is  removed 
in  a magnetic  field  and  that  the  elastic  anisotropy  of  crystalline 
compounds  can  be  nullified  and  even  inverted.  Most  of  the  observed 
effects  are  accounted  for  by  linear  magnetoelastic  theory  if 
b2>>t>i  and  proper  account  of  the  fluctuation  in  magnetization  is 
taken.  Discrepancies  indicate  the  presence  of  non  linear 
magnetoelastic  interactions. 
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